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' Abstract 

H 

Cd I Starting from an extension of the Poisson bracket structure and Kubo-Martin- 

Schwinger-property of classical statistical mechanics of continuous systems to spin 
systems, defined on a lattice, we derive a series of, as we think, new and interesting 
bounds on correlation functions for general lattice systems. Our method is expected 
to yield also useful results in Euclidean Field Theory. Furthermore the approach 
is applicable in situations where other techniques fail, e.g. in the study of phase 
transitions without breaking of a continuous symmetry like i-*((/))-theories with 
(hix) scalar. 
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1. Introduction 



In more recent times it has become more and more apparent that the natural Poisson 
bracket structure of systems of classical point particles is of considerable conceptual and 
practical value in the investigation of all sorts of questions being related to e.g. phase 
transitions in thermal equilibrium systems, transport equations and the like (see e.g. refs. 
[1] to [8]). 

One reason of its usefulness derives from the fact that both the classical equilibrium 
condition (more specifically, the so-called Kubo-Martin-Schwinger (KMS)-condition; cf. 
e.g. [9]) and the concept of symmetry (breaking) can be very neatly implemented by this 
structure. 

The KMS-condition in the regime of classical statistical mechanics of point particles 

reads: 

<{A,B}>^P-<B-{A,H}> (1.1) 

with A, B local observables (without loss of generahty: real, differentiable functions with 
compact spatial support on phase space), H the Hamiltonian, (3 inverse temperature, 
< • > the thermal average and the Poisson bracket given by 

{A, B} = J2idrA . d,^B - d,A ■ dr^B) (1.2) 

i 

(1.1) was already derived by Mermin (cf. ref. [10]) and was used by him to get the 
classical counterpart of the so-called Bogoliubov inequality 

< {A,B} >^ < p < A^ >< {B,{B,H}} > (1.3) 

which underlies, in various disguises, some of the approaches developed e.g. in the above 

cited papers. 

On the other side there exists a considerable amount of model systems in classical sta- 
tistical mechanics which do not openly carry such a nice structure, to mention typical 
cases in point, systems living on a discrete base-manifold as e.g. spin systems or systems 
defined via functional integrals. 

In view of the great calculational advantages of the Poisson bracket formahsm in classical 
statistical mechanics it would be highly desirable to have an analogous machinery at ones 
disposal for such (non-canonical) systems. 

The necessary general steps in this direction have been undertaken by us in ref. [11]. 
Furthermore we studied a (however rather limited) class of models (e.g. interface Hamil- 
tonian) to demonstrate the usefulness of our approach but made no systematic analysis of 
more complicated classes of spin systems which, while being of genuine relevance in sta- 
tistical mechanics as such, serve furthermore as the starting point of Euclidean quantum 
field theory. 

That is, we will derive in the following in a more systematic manner various (possi- 
bly new) classes of correlation (in) equalities for spin systems and apply them to typical 
problems in this field. 
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2. The Basic Estimates 



In this chapter we will apply our general results ([11]) to spin systems. Furthermore, to 
keep matters transparent, we will treat, in a first step, continuous spin systems, i.e. each 
spin S{x), X G Zi'^, ranges from plus to minus infinity. 

This is a sufficiently large class of models being frequently discussed in the literature. 
Our formalism can of course also be applied to spins living on, say, a compact manifold 
as e.g. S'^ etc. 

In that particular case the independent variables defining the configuration space would 
be certain angles. With {S{x)} — R the technical manipulations are a little bit simpler 
(for a discussion of e.g. the x — y— model along these lines cf. [11]). 

2.1 Remcirks: i) One can, at least in principle, get corresponding results for constrained 
spin systems as limiting cases by approximating their support, e.g. given by 5(5*^ — 1) • 
dS, S G R'^, with the help of a sequence of smeared out distributions, taken from som 
function space, Pn{S^ — 1), which converge in this limit towards 6{S^ — 1). 

ii) Technically it turns out to be advantageous to absorb a possible extra weight function 
p{S) occurring in the single-spin distribution into the Hamiltonian by writing it as an 
exponential. 

With the base space being Z'' or a • Z*^, a, the lattice spacing, statistical mechanics on a 
subset A C Z*^ is defined in the usual way via a Hamiltonian H{S_),S_ :— {S{xi)}, Xi G A, 
i.e.: 

<A> := Z-^ ■ j A{S)e~^"^^dS 
with dS_ := n^. ^ j\^dS{xi), 

Z := je-^^'^^MS (2.1) 

(possible extra weight factors being absorbed in H{S_)\). 

As to the class of admissible observables we choose them to be real, differentiable func- 
tions with respect to the variables S{xi) and, if necessary, bounded away from the internal 
and (or) external system boundaries in order to avoid artificial boundary terms in the var- 
ious partial integrations. This is however only a measure of precaution since, typically, 
at the internal system boundaries they are strongly supressed in most cases by the expo- 
nential vanishing of exp{—(3H) for some S{x) approaching ±oo. 

Now, as symmetries induced by Poisson brackets are locally generated by certain first 
order partial differential operators, we generalize the Poisson bracket structure of clas- 
sical point mechanics by replacing Poisson brackets with general first order differential 
operators acting on the spins in the following way: 

2.2 Definition: A first order differential operator operating on the spin system is given 

by 

V:=Y.d,{S,)ds, (2.2) 
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with di certain twice differentiable functions of the spin variable Si {Si abbreviation for 
S{xi)), denoting partial differentiation with respect to Si. V is acting on observables, 
introduced above, i.e. on certain differentiable functions of the spins {Si}. 

2.3 Remark: Note that the diS can in principle carry an extra dependence on the site 

i, where they are localized. 

With the help of such D's we are able to generalize the KMS-condition (1.1). 

2.4 Generalized KMS-Condition: 

<VA> = <A- ((3VH - div d) > 
withdivd := Y^dsM^i) (2.3) 

i 

and after some calculations the Bogoliubov inequality (1.2). 

2.5 Generalized Bogoliubov Inequality: 

<VA>^ << A>'^ ■ < (3- VVH - V{ div d) > (2.4) 

2.6 Remarks: i) That these are the proper extensions of the relations (1.1), (1-2) to 
lattice systems has been shown in [11]. 

ii) We have given another, slightly different extension in [11] via augmenting the local 
phase space at each site Xi and the Hamiltonian H by means of which we get a true 
Poisson bracket structure also on the lattice. 

Similar formulas can be derived for spin systems defined over a continuous base space, i.e. 
Euclidean quantum field theory over R'^, where expectation values like (2.1) are replaced 
by functional integrals and the derivative operator D by a functional derivative: 

< A >:= Z-^ ■ J L>[0]A[0]e-^^['^] (2.5) 

D[0] being the functional measure, H[(f)] the Euclidean action 

P := y" d'^x d{x, (j){x)) ■ 6/6(j){x) (2.6) 

(cf. [11]). 

As one has to cope in this continuum situation with various renormalisation problems 
if one does not treat the expressions in a purely formal manner we plan to study this 
situation in a more systematic way in forthcoming work. On the other hand one can try 
to carry over the corresponding expressions from the lattice situation by taking the lattice 
spacing to zero. 

2.7 Remark: In Euclidean field theory a certain method has been in use which bears a 
weak resemblance to our approach and which is called "The Integration by Parts Method" . 
(cf. e.g. [12]). Our method however appears to be considerably more general and yields 
stronger results since it draws on concepts which have not been exploited up to now in 
that field. 
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3. Correlation Inequalities for the (/)^-System 



The typical regime of application of relations of the above type is the important field 
of phase transitions and spontaneous symmetry breaking (examples can be found in the 
above mentioned literature). 

In that situation the derivative operator T) is typically chosen to be the generator of 
the flow representing the continuous "formal" symmetry of the model. The attribute 
"formal" means that in case the symmetry is spontaneously broken the Hamiltonian H 
is only invariant under the symmetry in a restricted (formal) sense and the manipulation 
of various limiting procedures becomes a highly delicate matter. 

In the following we want to show that, perhaps a little bit surprisingly, the above formulas 
can be applied also in a much more general environment, e.g. where no continuous 

symmetry exists at all or where the symmetry is unbroken. 

In this wider context the differential operator T) does not suggest itself but has to be 
chosen cleverly in order to yield interesting (in) equalities between various expectation 
values or correlation functions. 

3.1 Remcirk: In the following we will choose the differential operator V in such a way 
that div X> = holds in formulas (2.3), (2.4). 

To begin with we take as model Hamiltonian the 0^-Hamiltonian [J^ = 0) : 

= 1/2 X: JiASj + Y.^Sl + \Y,St (3.1) 

ij i i 

and discuss this model for various choices of the parameters. 

As interesting and more generic phenomena arise in the thermodynamic limit A — > Z'^ 
the following estimates are understood in this limit, i.e. they are calculated for finite A 
and are then generalized to the infinite system by a standard procedure. 

The different classes of estimates arise from different choices of the observable A and 
operator V in formulas (2.3), (2.4). The limiting equihbrium state is assumed to represent 
a "pure phase". This can be achieved in the usual way by adding a symmetry breaking 
term in H which is switched off in the end or by fixing certain boundary conditions. As 
both A and V arc chosen to be local, i.e. among other things, supported away from 
the system boundaries, there remains no explicit effect of these boundary conditions "at 
infinity" . 

3.2 The Case A -.^ Si,V := % 
Inequality (2.4) yields the estimate 

l^<l><P<Sf> ■{2m + 12A- < Sf >) (3.2) 

with the Hamiltonian (3.1). 
Solving for < Sf > we get: 
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i) m > : 



ii) m = : 



+ 7»2 _ ^ 



iii) m > 0, A — )• : 



< > > 77^ (3.5) 



2m/3 



iv) m < (i.e. ground state degenerated): 
Prom (3.2) we infer: 



r'<12A.«S?>+3^f-^ (3.6) 
and 

2m + 12A < 5^ > > always (3.7) 

which yields 

in 

< 5"^ > > -— always (3.8) 

Inserting this in (3.6) we see that the bracket is always positive also for m < 0!, i.e. there 
is no problem with the squareroot and we get, as in the case m > 0: 



<S!>> ^ ^ , m < (3.9) 

3.3 Remark: Note that our estimates hold both for the case < Si >— and < Si >^ 
0, i.e. with or without spontaneous magnetization. That is, < Sf > is either a pure 
fluctuation term or contains the overall magnetization < Si >. 

As an application of our above estimates we will study the case m < more closely. 
This is the regime where spontaneous symmetry breaking becomes possible. It is then 
an important question for which values of the phase space parameters, e.g. {P, m, A} the 
equilibrium state of the system represents the non-degenerated phase (i.e. single phase 
region) and for what values the equilibrium state is degenerated (i.e. two pure phases). 

3.4 Remcirk: At the moment the relation of our results (presented below) to other kinds 
of estimates of these bounds (sec e.g. [12]) is not entirely clear to us since they are usually 
derived by completely different methods. This point shall be clarified in the future. 

We now study the particular model Hamiltonian 

H = Wj-(Si-S,)' + Y.{moS^ + XSt), J>0 

nn i 

= -j2J-SiSk + J2i^s^ + ^st) (3.10) 
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with m :— rrio + n ■ J,n — number of nearest neighbors. 

For mo < 0, J > 0, i.e. ferromagnetic couphng, the ground state {(3 — oo), has S — 
±Y^| mo I /2A. With the help of our estimate for < Sf > we can provide a lower bound 
of the strength of fluctuation of Si and can set it into relation to the inverse temperature 
/9. That is: 

<S^> -S'^ > S'^ imphes 
<{Si-<Si>f> > s^=\j^ for < |< 5i >| <| :S I . (3.11) 

We can now estimate the critical value for (3 so that the mean deviation from the average 
< Si > becomes larger than |< 5'j >| or | 5" | itself. We conjecture that this signals the 
transition from the two phase to the one phase regime as typical fluctuations will then 
connect the two minima. (We, however, do not intend to prove this at this place). 

Prom (3.9) we get: 

3.5 Observation: (3.11) holds if 



^ TTTl o„ T I I ( and if /3* > 0!) (3.12) 



A 

14mo — 2nJ- \ mo 

i.e. we suppose that for P < j3* the system is in the one-phase regime. 
3.6 The Case A := ■ Sj, V = ds^ 
Inserting these expressions into (2.4) yields: 

< Si >'< P<S^-S]> (2m + 12A < S] >) (3.13) 

If the system is in a pure state one can exploit well-known cluster theorems (cf. e.g. 
[13]) to infer: 

< S^ ■ S] > — ^< S^ > ■ < S] > ^ < Sf >^ for I i - j H oo (3.14) 

This yields: 

< >2< (3{2m < Sf >2 +12A < Sf >'^) (3.15) 



and for m < 0: 



3.7 The Case A := Sf ■ Sj, V = dg^ 

< Sf >'< P-<Sf- S] > •(2m + 12A < S] >) (3.17) 
and with the cluster property: 



It is evident that one can derive, proceeding in the indicated manner, a whole sequence 
of inequahties between various expectation values and correlation functions. 

3.8 Remarks: i) Up to now our estimates are in general independent of the strength 
of the ("kinetic") coupling Jij. This is a consequence of the choice S^. for V. For T) — 
dsi + ds^.i 7^ j- WH yields also terms containing the couplings Jij. 

ii) Furthermore, our results are dimension independent. This is, however, not always 
a disadvantage. The dependence on the space dimension has, on the other side, been 
exploited in previous work of us (cited above). To incorporate dimension one has to 
choose observables A which go with the volume A. In this paper A was fixed independent 
of A. 

3.9 A Certain Strategy: 

Our estimates yield bounds from below, i.e., as a case in point: 

< 5"^ > > expression (1) in (/?,m. A) (3.19) 

If it is possible to derive a bound of the sort: 

<{S- < S >)2 > < expression (2) in {l3,m, A) (3.20) 
it may become possible to get interesting bounds of the kind: 

< 5" > 7^ for a certain regime of (/3, m. A) (3-21) 
i.e. estimates concerning the existence of phase transitions. 

3.10 Observation: Note that estimate (3.8) is even independent of the inverse temper- 
ature /3, i.e. 

<Sf>> ^—^ for m < (3.22) 
' 6A 

This implies that in the limit /3 — > oo we end up in a completely ordered phase, i.e. spins 
aligned, or, in the cases where an ordered phase is excluded by e.g. Mermin- Wagner- 
theorem, the ground state consits of a random occupation of the two minima of the 
Hamiltonian. 

In closing this paper we would like to point to the fact that up to now we have only 
exploited the inequality (2.4). Inserting, on the other side, the various choices for A,!) 
into (2.3) we get an equality between certain expectation values. For the simplest choice, 
i.e.: 

3.11 The Case A := Si,V -.^ ds^, Equation (2.3) 

f3-^ =< J2 JijSiSj + 2mSf + 4XSf > (3.23) 
j 

which is reminiscent of sort of a virial theorem for spin systems. Corresponding equations 
can be derived for other choices of A, V. 

3.12 Remark: While we think the above scheme is representing a new approach to this 
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field it may well be that the results can possibly also be deduced by employing other 
methods. In any case, an advantage of our approach is, in our view, its simplicity and 
transparency. 
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